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SUMMARY 



A method for calculating the rate of heat transfer 
from the surface of an airfoil or streamline body is pre- 
sented, A comparison with the results of an experimental 
investigation indicates that the accuracy of the method is 
good • 

This method may "be used to calculate the heat supply 
necessary for heat de-icing or in ascertaining the heat 
loss from the fuselage of an aircraft operating at great 
altitude, for example. 

To illustrate the method, the total rate of heat 
transfer from an airfoil is calculated and compared with 
the experimental result. 



INTRODUCTION 



During the design of the test model for an investiga- 
tion roportod in reference 1, it was necessary to extend 
the theory of heat transfer into the laminar flow along a 
flat plato (rofcroncG 2) to cover the problem of determin- 
ing tho local r^^te of heat transfer from the surface of an 
airfoil into its laminar boundary layer. Using the same 
general method of attack, a solution was also obtained for 
tho turbulent flow region, making possible tho calculation 
of the total rate of transfer from the airfoil. 

Since tho boundary-layer characteristics of a stream- 
lino body may bo calculated in the samo manner as for two- 
dimensional flov7 if the spreading or crowding together of 
tho boundary layor duo to tho change in body dimensions 
along the axis is considered, tho heat-transfer equations 
developed for an airfoil were extended to make possible 
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tho computation of the total heat loss of a stroainline body. 
The derivation of tho method is given in tho appendix. 

A^limitod oxporimontal invostigation of the method 
was made. Those tests were conducted with the same test 
model as for reference 1 to determine the accuracy gf th:b 
method in calculating tho local rate of heat transfer into 
both the laminar and the turbulent boundary layer of an 
airfoi.lt and also to obtain a check on the computed total 
rate of heat transfer for the win.?;. 

It is hoped that this method will facilitate a more 
accurate determination of the heat losses from wings in 
designing heat de-icing systems as well as fr©m fuselages 
in the design of cabin-heating systems for aircraft oper- 
ating at great altitudes. 



SYMBOLS 



0?he symbols used throughout this report and in the 
appendix nre defined as follows: 

V free-stream velocity 

U local velocity just outside the boundary layer 

u local velocity inside the boundary layer 

c wing chord 

L length of streamline body 

X- distance along the chord from the leading edge for 

an airfoil, or along the axis for a streamline body 

y distance normal to the surface 

8 distance along the surface from the stagnation point 

r radius to surface of streamline body at any point 
along the axis 

p air density 

\X' absolute viscosity 



kinematic viscosity, ix/p 



surface shear 




tur"buleiit toundary-layer parameter, ^ = 



Reynolds number based on wing chord, Vc/u 

Seynolds number based on body length, VL/u 

local Mach number, rntio of velocity just outside the 
boundary layer to local velocity of sound 

free-stream Mach number, ratio of velocity of free 
stream to the velocity of sound in the free stream 

momentum thickness of boundary layer, 



heat-transfer characteristic length for a laminar 
boundary layer 

heat-transfer characteristic length for a turbulent 
boundary layer 



free-stream air temperature, rjbsolute 

local temperature outside the boundary layer, °F 

local temperature outside the boundary layer, °P abso- 
lute 

local temperature inside boundary layer, ^5" 

local temperature inside boundary . layer , absolute 

surface t emperatur e 

surface temperature corrected for compressibility 
-t^) heat-transfer temperature difference 

specific heat at constant pressure 



h 



e 




free-stream air temperature, 
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k heat conductivity 

y ratio of specific heats 

a Prandtl number, Cp|J./k 

e eddy viscosity 

P eddy heat conductivity 

heat -transfer coefficient, Btu/sq ft, ^F, sec 

q^. local rate of heat transfer, Btu/sq ft, sec 

Method 

The detailed ahalysis given in the appendix develops 

the following formulas, ty which tho local rate of heat 

transfer into hoth turbulent and laminar- boundary layers 

may "bo computed. This method. is applicable either to an 
airfoil or to a streamline ' body • The local rate for lam- 
inar flow is 



q^^ = 0.700 ^ (tp ^ t-o) (l)i 



or the local heat-transfer coefficient is 



hx = 0.700 (l)b 



and for turbulent flow, 



CLx = 0-760 ^'(tp - to) (2)a 



or the local heat-transfer coefficient is 



hx = 0.760 — - (2)b 



H eat transfer from an airfoil ,- For the laminar 
boundary layer of an airfoil, 8^ . is- computed as in ref- 
erence 3 from the pressure distribution, as follows: 
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=L = c / ^ TS^V / ^ d f (3) 



For the turbulent "boundary layer of an airfoil i 
is computed as 

6^ = — (4) 



where I is the value of the turhulent toundary-layer 
parameter as determined l^y a st ep-"by--step solution of the 
relationship of reference 4, given as 

dl 6a3dU^U 

.dx ^ U ix ^ ^p ^^^^ ^""^ 

The value of f ( ^) is given in table I as taken from 
refcronco 4, and may "be plotted on semilogar i thmic paper 
for ease in using. 

Eeat transfer from p. streamline body For the lami- 
nar boundary layer of a streamline body, 6^ is computed 
from reference 3 as 



The turbulent boundary-layer heat-transfer length may 
be computed as 




8a. = —Lk-. (7) 



where ^ is determined as for an airfoil from the equa- 
tion 

dt ^ 6.13 dU , 2.557 dr U ./ /^v 
di * d^ ^ di = V ^^^^ 
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"by the st ep-lDy- st ep process mentioned for calcvilating Sq; 
for an airfoils 

C om-pres silplllty c or rect i o n If the heat flow is to "be 
obtained at free-stream Mach numbers such that the aerody- 
namic temperature rise is an appreciable portion of the to- 
tal temperature difference, a correction for aerodynamic 
heating should he made. The "heat-transfer temperature 
difference" to "be used for a laminar boundary layer is 

(tp - to) = (t, - t^) - 0.20 M^^ fl -^0.13 (^JJ (9) 

and for turbulent flow 

(tp - t^) = (t^ - t^) - 0.20 To (10) 

where (t^ - t^) is the desired temperature rise for 
heat de-icing. 

The total rate of heat transfer from an airfoil or 
streamline body may be found as follows: 

1. Estimate the location of the transition point by 

the method of reference 5 for an increasing 
pressure gradient, or by reference 5 for a 
falling pressure gradient. 

2. Calculations for the laminar region ahead of the 

transition point . 

a) Compute the values of 6^ along the surface 

to the transition point by equation (3) for 
an airfoil, or by equation (s) for a stream- 
line body. 

b) With those values and the desired temperature 

distribution corrected for compressibility, 
compute the local rates of heat transfer 
along the sixrface by equation (l)a. 

3. Oalc^ilat ions for the turbulont region behind the 

transit ion point . 

a) S'rom the value of 6l at the transition point 
compute 9» the -momentum thickness as 



e = 0.289 6l 
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Using this vrIug of 6, find tho initial 
valuG of I at the transition point as 

i = 2.557 logQ (4.057 

Id) With this initial value of ^, calculate the 
values of along the surface by equation 

(5) for an airfoil, or by equation (8) for 
a streamline body* With these values of , 
compute 851 along the surface by equation 
(4) for an airfoil, or by equation (7) for 
a streamline bod;/. 

c) Using these values of 8^ -'^nd the desired 

t emper .c? ture difference across the boundary 
l^yer corrected for compressibility, com- 
pute the local rate of he^t transfer along 
the surface by equation (2) a. 

4, Integrate theso local rates 0|f heat transfer along 
the chord for both laminar and turbulent regions 
to obtain the t otal rate of heat transfer. 



Heat -Transfer Measurement s 

As mentioned in the Introduction, a limited number of 
heat-transfer tests were made on the heated wing model of 
reference 1 in the 7- by 10-foot wind tunnel of the Ames 
Aeronaut icp.l La.bor-itory to check the accuracy of the theo- 
retical iiiothod. 

Since the wing was not designed specifically for heat- 
tr'^.nsfcr tests, the experimental results r.ro subject to 
sGvor.al sources of error. All computed ratoe 01 heat 
transfer are based on the t ei-aporaturc distributions ob- 
tained at tho center of tho 3pan» -assuming that tho span- 
wise variation is nogligiblo. This ic* escontially true 
except for a small portion at each end of the wing. Pre- 
cautions were b^ikon to minimize tho heat losses at tho 
ends of the wing. Those r.TQ not boliovod large since tho 
design of the heating system allowed only slight transfer 
by convection, nnd the conduction of heat from the wing 
to its supports is negligible. Losses due to radiation 
from the v/ing have been computed as a maximum of 5 percent 
fsr the whole surface heated to 100° F above the surround?- . 
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ings. This loss is not considered in the he?>.t-t ransf er 

d^^ta^ 

For the purpose of computing heat transfer, the chord- 
v;ise temperature distribution v/as computed ^nd plotted as 
"heat-transfer temperature difference," (tp - t^). This 
is the observed temperature difference corrected for com- 
pressibility effect. That, is, tp is the temperature 

measured by a thermocouple in the skin minus the computed 
aerodynamic heating temperaturo rise. The value t^ is 
the free-»strcam air temperature (distinguished from the 
local temperature just outside the boundary layer of the 
v/ing, x^rhich \irill be higher or lower than t^ due to adia- 
batic variations caused by the velocity field of the wing). 

RESULTS AHD DISCUSSION 



Heat-transfer tests of the wing were made in two 
parts, the first concerned \\^ith comparing computed and 
measured values of the local rnte of transfer from the air- 
foil surface to laminar and turbulent boundary layers, 
and the second with checking the total rate of heat trans- 
fer from the airfoil. 

The tests to measure the local rate of heat transfer 
were i^cade at Cj = 0 for two test Reynolds numbers, with 

free transition to obtain the heat-flow rate into a lami- 
nar boundary layer, nnd with transition fixed at 5-percent 
chord to determine the flow rate into a turbulent boundary 
layor» Figure 1 presents the pressure distribution over 
the v/ing at Cj = 0. The " experiment al procedure consisted 
in adjusting the heat input so that the skin temporaturos 
\TOTQ nearly constant nlong the chord. With this tempera- 
ture distribution achieved, it was assumed that the second 
compartiaent of the i^ing, extending from 14.6-percent to 
26.3-percent chord, was thermally isolated so that no* flow 
of heat occurred in the skin or through the bulkheads. 
The power input to this comp-^rtment wa? then measured by 
raeans of a voltmeter and an arameter for comparison with 
the calculated rate of heat flow. 

Figure 2 shows that the desired constant chordwise 
temperature distribution v;as attained for the laminar 
boundary layer, but the data of figure 5 for . turbulent 
flow indicate that while the distribution was nearly con- 



9 



strmt from 10- to 30-percent chord, covering the region 
under consideration, the temperatures over the nose were 
excessively high. This came about through the heating dif^. 
ficultios resulting from the sudden change in heat-transfer 
coefficient at the point where transition was fixed. This 
typo of distribution may have resulted in some change in 
the local values of the temperature gradient at the wing 
surface, though this effect should be small since no ap- 
preciable temperature gradient existed over t}ie portion of 
the surface concerned. 

To obtain the computed values, the variation of the 
heat-transfer coefficient along the chord was calculated 
for each case as outlined under the section Method (p. 4) 
(results plotted in figs. 4 and 5). The heat input into 
the second compartment was then computed from the heat- 
transfer coefficient and experimentally measured tempera- 
ture difference (corrected for compressibility heating 
effects)* Both the measured and the computed values of 
heat input into the second compartment are listed in the 
follov/ing table: 



a 
(deg) 




Boundary- 
layer 


Calculated 
heat trans- 
fer, (kw) 


Measured 
input , 
(kw) 


Err or 

input , 
( percent ) 


0 


6.71 


laminar 


0.861 


0.875 


1.6 


0 


10.72 


■_-do--- 


1.048 


1.038 


1.8 


0 


6.90 


Turbulent 


3.76 


3 .75 


.3 


0 


11.17 


__do 


4,01 


4.01 


.0 



The check is considered quite satisfactory, and is taken 
to indicate that the method for the computations of the 
heat-transfer coefficient involves no serious errors de- 
spite the assumptions involved. 



Further tests for the purpose of establishing the 
validitjr of the method as regards the total rate of heat 
flov/ from a wing were made at a lift coefficient of 0.55 
and 8.60 million Reynolds number. In addition to a test 
with free transition, a second condition simulating the 
formation of ice near the stagnation point by fixing tran- 
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sition at 5-percent chord on the lower surface was invest!- 
gntGd» Hoat input in "both conditions was maintained at the 
maximum availa'ble from the appara,tus, and no r^.ttempt v/as 
made to achieve a predetermined temperature rise or chord- 
wise distributiono The chordwise temperature distribution 
obtained, corrected for the effects of compro ssibill ty 
heating, is plotted in figures 6 and 7. 

The chordwise variation of hoat - t ransf er coefficient 
was computed for both the transition-fixed and transition- 
free conditions, as outlined in Method. Tables II and III 
present the computations for the transition-free condition 
and serve as an illustrative ezample* The pressure distri-^ 
bution used for these calculations is given in figure 8. 
The value of heat input was then computed by use of the 
experimental temperature distributions (corrected for cam- 
pre s sibility) nnd the calculated heat-transfer coefficients 
of figures 9 and 10. The computed and measured heat input 
are compared in the table below. Results indicate satis- 
factory agreement. 



a 

(deg) 


R(,X10~®. 




Boundary 
layer 


Calculated 
heat 
transfer 


Measured 
input , 
(kw) 


Error 
input, 
(percent ) 


5 


8.69 


Upper surface 
turbulent 

Lower surface 
laminar 


19.60 


23.52 


16.8 


5 


8 .60 


Upper surface 
turbulent 

Lower surface 
turbulent 


22,50 


23 .55 


5.2 



The experimental temperature distributions for the 
tests of the total' rate of heat transfer show that the 
heat-transfer temperature difference varies to a marked 
dogroo along the chord (figs. 6 and 7). This variation 
violates one of the assumptions underlying the develop- 
ment of the method; that is, that the temperature differ- 
ence is constant along the chord, v/hich must bo true if 
tho transfer of heat at all points along the surface is 
analogous to the. transfer of momentum. To what degree 
this assumption may bo ignored has not been determined 
analytically since the problem of considering tho variation 
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of tcriiperaturo along tho chord presents difficulties which 
have no far prevented a solution. The expor imental re- 
sults, however, indicate that tho accuracy with which the 
total ro.to of heat transfer can he computed is not greatly 
impaired "by the temperature variations experienced. Oenor- 
alization of- this result must await further experimental 
checks, 

The accuracy with which the local rate of heat trans- 
fer may he computed in a falling pressure gradient is de- 
pendent upon the accuracy with which the surface shear may 
he determined. Squire and Young's method (reference 4) as- 
sumes that the turhulent houndary layer* in a falling pros- 
sure gradient exhibits tho same characteristics as the 
fully developed turhulent layer of a flat plate. The ex- 
tent to which the r olat i onship ho tween the surface shear, 
the momentum thickness, and the local velocity so derived ' 
remains valid is shown by the accuracy of the Squire and 
Young method in determining friction drag. It must he 
realized, however, that tho method will fail if turhulent 
separation is imminent.- 

The thickening of the turbulent boundary layer due to 
the fall, in pressure acting on the displaced mass of fluid 
pIso is ignored by the assumption that tho heat -transfer 
rate is proportional to the surface shear computed by 
Squire and Young's method. Actually, the heat capacity of 
the boundary layer is incre?^sed by this thickening which 
tends to increase the rate of heat flow at the surface. 
This counteracts tho effect of the profile distortion, re- 
sulting from the same cause, v;hich roducos the surface 
shear since it tends to cause separation. But this ef- 
fect, too,. is negligible for all cases where Squire and 
Young's method may be applied. 

In concluding, it must bo stated that while the meth- 
od presented herein is subject to a number of broad as- 
sumptions in its development, tho experimental evidence 
presented shows the total rate of heat flov; may be calcu-. 
lated v/ith reasonable accuracy. 



CONCLUSIOFS 



The accuracy of tho method for determining the rate 
of hea-t transfer from an airfoil is shown t o be good by 
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tho results of a limitod oxper imontal investigation. 
Since the corroctnoss with which tho heat transfer can ho 
computed is dependent mainly cn the accuracy with which 
tho boundary-layor characteristics may ho detorminod, it 
is expected that the method possesses the same accuracy 
for computing hoat-t ransf or rates from a streamline "body. 

Although tho development of the heat -transfer formu- 
las is based on the assumption th^t tho skin temperature 
remains constant along tho surface, the experimental re- 
sults show that for moderate temperature variations the 
precision is still good. 



Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 



1 • Hea.t transfer into a laminar boun d ary layer ; - The 
theory of heat ttransfer into a laminar boundary layer was 
first investigated by E. Pohlhausen for the case of incom- 
pressible flow along a flat plate maintained at a constant 
temperature (reference 2). Pohlhausen* s solution is de- 
veloped by solving the differential equation for the tem- 
perature boundary layer by using Blasius' solution for 
the velocity boundary layer. 

In order to arrive at a solution for an airfoil in an 
incom-pre ssible fluid, it is necessary 

(l) to assume that the temperature boundary-layer and 
the velocity boundary-layer profiles for the airfoil are 
related in the same manner as for Pohlhausen' s solution. 
(This is true if the temperature of the skin remains con- 
stant along the surface and -.if the thinning of the fric- 
tion layer in a favorabl e pressure gradient due to the 
change in pressure acting on the displaced mass of fluid 
is negligible . ) 



APPENDIX 



locity boundary layer and then determine 



(2) to calculate the value of 
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of tcmperaturo along tho chord presents difficulties which 
have BO far prevented a solution, Ihe experimental re- 
sults, however, indicate that the accuracy with which the 
total rate of heat transfer can be computed is not greatly 
impaired by the temperature variations experienced. &ener- 
alization of this result must await further experimental 
checks. 

The accuracy with which the local rate of heat tra'ns- 
fer may be computed in a falling pressure gradient is de- 
pendent upon the accuracy with which the surface shear may 
be determined. Squire and Young's method (reference 4) as- 
sumes that the turbulent boundary layer in a falling pros- 
sure gradient exhibits the same character i st i cs as the 
fully developed turbulent layor of a flat plate. The ex- 
tent to which the relationship between the surface shear, 
the momentum thickness, and the local velocity so derived ' 
remains valid is shown by the accuracy of the Squire and 
Young method in determining friction drag. It must bo 
realized, however, that the method will fail if turbulent 
separation is imminent. 

The thickening of the turbulent boundary layer du'e to 
the fall- in pressure acting on the displaced mass of fluid 
p.lso is ignored by the assumption that the heat -transfer 
rate is proportional to the surface shear computed by 
Squire and Young's method. Actually, the heat capacity of 
the boundary layer is increased by this thickening which 
tends to increase the rate of heat flow at the surface. 
This counteracts the effect of the profile distortion, re- 
sulting from the same cause, v;hich reduces the surface 
shear since it tends to cause separation. But this ef- 
fect, too, is negligible for all cases where Squire and 
Young's method may be applied. 

In concluding, it must bo stated that while the meth- 
od presented herein is subject to a nuraber of broad as- 
sumptions in its development, the experimental evidence 
presented shows the total rate of heat flow may be calcu- 
lated with reasonable accuracy. 



CONCLUSIOirS 



The accuracy of the method for determining the rate 
of heat transfer from an airfoil is shown t o bo good by 
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tho results of a limited oxpor imont al investigation. 
Since the correctness with which the heat transfer -can be 
computed is dependent mainly on the accuracy with which 
tho boundary-layer characteristics may bo detorrained, it 
is expected that the method possesses tho same accuracy 
for computing heat-transfer rates from a streamline body. 

Although the development of the heat-transfer formu- 
las is based on the assumption that the skin temperature 
remains constant al ong the surface, the experimental re- 
sults show that for moderate temperature variations the 
precision, is still good. 



Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Koffett Field, Calif. 



APPSNDIX 



1« Heat t ransfer into a laminar boundary layer ,- The 
theory of heat ttransfer into a laminar boundary layer was 
first investigated by E. Pohlhausen for the case of incom- 
pressible flow along a flat, plate maintained at a constant 
temperature (reference 2). Pohlhausen' s solution is de- 
veloped by solving the differential equation for the tem- 
perature boundary lay.er by using Blasius' solution for 
the velocity boundary layer . 

In order to arrive at a solution for an airfoil in an 
incom-pressible fluid, it is necessary 

(l) to assume that the temperature boundary-layer and 
the velocity boundary-layer profiles for the airfoil are 
related in the same manner as for Pohlhausen* s solution. 
(This is true if the temperature of the skin remains con- 
stant along the surface and if the. thinning of the fric- 
tion layer in a favorable pressure gradient due to the 
change in pressure acting on the displaced mass of fluid 
is negligible.) 

(3) to calculate the value of for the ve- 

locity toundftry layor and then doterinine ( , ) with 
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the relationship resulting from (l). (The solution of tho 
prohlom for the temperature of tho skin varying along tho 

chord has "boon prevontod because tho difficulties so far 
have hocn found insurmountable.) 

Pohlhausen*s expression for the temperature gradient 

in the "boundary layer at the- surface of the plate is given 

as 



The function a(a) is the first derivative of Pohl- 
hausen's function defining the temperature boundary layer 
which, for a = 1, is equivalent to the second deriva- 
tive of Blasius' function for the velocity boundary layer, 
Pohlhausen found, that a(or) is accurately given by the 
relationship 





then 




Now, for 



the Blasius boundary-layer distribution 




so that 




(If) '-/^ 

y=o 



(^p - ^o) 



y=o 



It is now necessary to determine 




for the airfoil 



at p.ny chordwise position. For the Blasius profile 




where 8l is the thickness oT the boundary layer where 
u = 0*707 U. Substituting in equation (2) 
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(~) = 0.765 /F — ~ 

or, taking a = 0.760 for air, the local rate of heat trans 
fer is 



a 



= ^ (-di),.^ = 0.700 A (tp - t,) (1) 

or the heat-transfer coefficient is 

fix « 0.700 ^ (l)t 

The development in reference 6 of an expression for 
the heat-transfer rate "based on Eeynolds analogy gives re- 
sults which are in complete agreement with the ahove if 
0=1. However, the experimental results of reference 7 

indicate that the expression ^0.760, that is, *J ^^ix* 

properly relates the velocity and temperature gradients in 
the laminar layer. The two methods give results within 10 
percent of each other, which is sufficient for practical 
cases. 

The values of 6^ for laminar flow may he determined 
"both for an airfoil and a streamline "body by the method of 
reference 3. 

2. Turhulent boundar y layer The theory of heat 

transfer in eddying flow as given by Dryden (reference 8) 
requires the introduction of several nev/ concepts. If the 
equations of motion for turbulent flov; are x^ritten by plac- 
ing u = u+u*, v = v+ v*, w = w + \^r', where the bars 
indicate mean values and the primes indicate fluctuations, 
and these values substituted into the equations of motion 
for steady flow, similarity between equations so developed 
and the stoady-flow equations can bo shown by introducing 
a value of eddy viscosity, €. Similarly, the concept of 
eddy^heat conductivity, P, is introduced by placing 

t = t + t' in the equations of the temperature field. 

Those values of eddy viscosity, c, and eddy conduc- 
tivity, P, ho\i^evor, do not have the same properties as 
and k since they vary from point to point in the flow. 
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Nevertheless, e and P can "be shown to vary in tho same 
manner from point to point in the fluid* This is done hy 
introducing Prandtl's concept of a mixing length; that is, 
a length of path follov/ed by a fluid particle before it 
becomes lost in the mass of eddying fluid. 



It is therefore shown if the shear, 



r 



then 



s= -pu * V * , 



or that 



pi 



du 


du 


dy 








dy 





in which the mixing 
in the fluid. Nov/, 



ered to be 



-CpV 



length, it varies from point to point 
if the eddy heat transfer is consid- 

then the eddy heat conductivity is 



ecLual to c-.pr 



or equal to 



Cp€, 



provided the mixing 



as.. 



du 
dy 

length for heat transfer is considered to be the same 
the mixing length for the transfer of shearing stress. 
Dryden states that available experimental data show that 
the mixing lengths near a wall are closely equivalent for 
transfer of heat and momentum, but that the relationship 
falls down, for instance, in the wake of a heated body. 
Since the present case concerns heat transfer from a wall 
to eddying flow in a boundary layer, it is believed that 
this relationship is acceptable. 



For turbulent flow, it has boen shown that Jhe Prandtl 
number is equal to unity; that is, 



CpM. 



turb _ 



Hurb 



If the Prandtl number is unity, then the thermal and dynam- 
ic boundary layers havo the same profile (reference ?)• If 
we make the same assumptions as in step (l) for the laminar 
boundary layer, v/o may write 
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where U is the velocity outside the friction layer. (This 
relationship is dependent on tho assumption that the tem- 
perature along the surface remains constant as for step 
(l) for the laminar boundary layer,, and that the thickening 
of the "boundary layer due to. the increasing pressures act- 
ing on the displaced mass of fluid and the distortion of 
the profile thus resulting is negligible.) 



but P ^p € 

so that = -cp€ (^) 



(tp - to) 



dyj U 
y=o 



since the surface shear 



Vdv^ 



^^y'y=o 



= Cp To 



U 



This is the same formula developed by fieynolds for flow 
in pipes (reference 8). 

So it is seen that the problem of calculating the rate 
of heat transfer in turbulent flow is primarily a prcblem 
of calculating the surface shear along the airfoil. This 
may be done by the method of reference 4, in which Squire 
and Young write the relationship 

where t,^ = 2.557 logg ^4.075 ^) 

8 "being the momentum thickness sf the 'bcundary layer. 
(This relationship is developed from van Karman's formula 
for tho skin "frict ion experiencod hy a flat plate with a 
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fully dovcloped turbulent boundary layer. This assumption 
"becomes loss and less true as the turbulent boundary-layor 
profile of an airfoil becomes' changed in shape and ap- 
proaches separation in a steep pressure recovery.) Substi- 
tuting for Tq 

pU^ ( t^ - t,) 

<lx - Cp -JS- — ^ u . 

pU . 

'Ix = ^'p "p" ^*p - 

or 

k U p, V c . N 

•^x " °pP k ^2 ]I V c " 

CpM- k Vc U - , , 
k V V P 

whoro M. and k aro valuos for unoddying flow. Sutsti- 
CpM. 

tuting for the value for air, 0.760, 

= 0.760 Be (j) (tp - t,) 

or considering — = 8t • a characteristic length for 



the turbulent boundary layer, 



then , = 0.760 -^ (tp - to) (2)a 

and hx = 0.760 (3)b 



6rp 

In calculating 6ij , ^ may be computed by the step- 
by-step solution of the equations of reference 4 for an 
airfoil or streamline body. 

3. Compressibility effects on heat transfer . Since 
the foregoing analysis has been made for incompressible 
flow, the effect of aerodynamic heating must be dealt with 
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if the heat transfer is to be accurately obtained. The ef- 
fect of compressibility may be considered simply as influ- 
encing the heat-transfer temperature difference to be used 
in the above-developed equations; that is, a part of any de- 
sired increase in the skin temperature will result from 
aerodynamic heating, and this part of the temperature In- 
crease involves no expenditure of heat, 

0?ho temperature field near the heated surface of an 
airfoil or body of revolution operating at high Mach num- 
bers may bo determined by superposing the heat-transfer 
temperature field on that due to the friction heating as 
in reference 9, Sckert (roforenco lO) has shown that the 
temperature field duo to aerodynamic heating for a = 1 
may bo o.ifprcssed as 



t = tL + 



V - 1 s-,, 



1 - 



'u 




or 



Superposing this on the temperature field for heat 
transfer, which may bo given as 



tho comljinod tomporaturo fiold is 

''UN -V _ 1 2, 



t = t 



V 



^ "i^^'^ Iff 



It is evident that at y = 0 



\dy/ 



(t. 



which indicates that the heat transfer corresponds to the 
heat-transfer temperature fiold, so that for compressible 
flow the only correction no-ccssary is thpt of correcting 
the skin temperature for tho rise in temperature due to 
aerodynamic heating. 



Eckort has shown that 



^-0^^ M/.,[x-(x -y.)(|;] 

We may write for the laminar region, <j = 0,760 

and for the turbulent region, a = 1 

(tp - to) = (t^'-.to) - 0.2 m/To 
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TABLS I 

fiO = 10.411 t 



Numerical values 



i 


10® f(^) 


i 


10^ f(t) 


14 


221.6 


24 


1.502 


15 


130.39 


25 


.937 


16 


77.47 


26 


.585 


17 


46.11 


27 


.367 


18 


27.97 


28 


.2308 


19 


16.96 


29 


.1454 


20 


10.35 


30 


.0917 


21 


5.35 


31 


.0582 


22 


3.906 


32 


.0369 


23 


2.418 


33 


.0235 
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TABLE II 

COMPUTATION OP HEAT-TEANSPBR COEFFICIENT 
Airfoil 65,2-016 c-j=0.55 £(.=8.69x10® ' Mo=0,lg3 



Upper surface Transition at chcrd 615^0=15^.5x10"® tl5fc=18.25 



Positien 


Ax 


(?) 


htt/v 
dx 


f(S) 

xxu 


d t 

CLX 


AS 




10^6-r 


1 

hx 


X 

c 


j 

1 (ff) 


0.01 


0.07 




1.61 










18.25 


166.5 


0.0191 




.lU 


0.07 




-0.667 


zk.e 


51.8 


3.63 








.02 




1.57 




1+.17 


10.36 




21,88 


21+5.1+ 


. .0130 






.07 


1.5U 


-.572 


. -73 








.03 


.21 












22.61 


267.1 


.0120 






.07 




-.41+3 


2,92 


7.3^ 


.51 






.28 




1.52 










23.12 


283.2 


.0113 






.07 




-.300 


2.29 


5.53 


.39 








.05 


.35 


,1k 


1.50 






I+.62 




23.51 


296.9 


.0108 




.^9 




-.270 


1.89 


.65 








.07 


.Ik 


1.1+8 










2l+.l6 


317.1+ 


.0100 




.63 


iM 


-.185 


1.39 


3.32 


.1+6 






.09 


,1k 










2I+.62 


335«3 


.00952 






1.I422 


-.133 


1.12 


2.57 


-36 








.11 


.77 












2l+,98 


353.^ 


.00903 


.16 




.35 




-.069 


,SkC 


1.96 


.69 






1.12 




1.399 










25.67 


379.5 


.0081+1 




l,kf 


'35 




-.01+6 


.685 


1.39 


.1+9 




.21 




1,382 




.5^ 






26.16 


399.0 


.00800 


.26 




.35 




-.0372 


1.10 


.39 






1.82 




1.386 




.1+53 






26,55 


1+10.0 


.0077s 






.35 




-.0286 


.906 


.32 






.31 


2.17 




1.356 








.26 


26.87 


1+30.9 


.OO7I+I 


.36 




.35 




-.0171 


.393 


.739 




1+39. 8 




2.52 




1.350 


-.O1I+3 


.3^5 


.61+3 




27.13 


.00725 


M 




.35 


1.3^ 


.23 






2. 87 












27.36 


1+1+8.2 


.00712 


.it6 




.35 




-.031+3 


.312 


.677 


.21+ 






3.22 




1.330 








27.60 


46l,0 


.00682 






.35 




-.0915 


.278 


.918 


.32 








.51 


3.57 














27.92 


1+72.1 


.00676 
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TABLE III 

COMPUTATION OF HEAT-TEANSPEE COEFFICIENT 
Airfoil 65,2-016 Ci=0.55 Ec=3. 69x10 Mo=0.1g3 

Lovxer surface Transition at minimum pressure 



X 


(i) 


6 ♦ 1 7 

(?) 


ci. / r* 

s 

'O 






10^ ii 

G 


10'6l 

(ft) 


10\ 


0.010 


0.0 






— — 


— 


— — 


— - 





.020 


.185 


o.ooogif 




: -i*xio~^ 


o-.oooi^-i 


0.775 


5.42 


•51^.2 


.030 


.290 


.0060 




3i+xio~® 


.0032 


.809 


5.66 


52„0 


.OlK) 


.370 


.0165 




lUUxio"^ 


.0100 


.9^+2 


6.59 


ii-9.6 


.050 


.^55 


.0388 




Ui6xl0r6 


.0261 


.992 


6.91+ 




.075 


.60 


.122 




23^0x10" s 


,09^ 


1.231I 




3^,1 


.100 


.71 


.21^3 




0.0056 


.205 


1.293 


9.06 


29.7 


.150 


.25 


.512 




.02U8 


.U72 


1.767 


12,37 


23.8 


.200 


.92 


.70s 




.0572 


.680 


2.267 


15,86 


ig<,5 


.250 


1.00 


1.00 




.1012 


1,00 . 


2,^4-95 


17,^6 


16.9 


.300 


1.07 


1.30 




.1610 


1.3'+ 


2.710 


18.97 


15.1 


.350 


1,12 


1.60 




,23^0 


1.69 


2.9U0 


20,i47 


1U.3 


.1400 


1.16 


1.81 




.3272 


1.95 


3.20 


22. IW) 


13.3 


.U5O 


1.17 


1.93 




.1^306 


2.10 


3.5^ 


2U.80 


11,8 


.500 


I.IS 


1.9s 




.5392 


2.15 


3.91 




11.0 . 
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Figure 5.- 'Chordwise distriljution of heat transfer coefficient, 
hx,a=0o, transition fixed at djl> chord. 
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Figure 6.- Chordwise distribution of temperature difference, (tp-.to)ia=50, 
transition free. 
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